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SECTION  1 
INTRODUCTION 


The  general  problem  to  be  considered  here  is  acoustic  scattering 
in  the  resonance  region  from  arbitrarily  shaped  3-dimensional  targets  in  a 
fluid  medium,  possibly  in  the  presence  of  a  bottom  (see  Fig. — ±}  .  The  tech¬ 
niques  presented  in  this  report  can  be  applied  to  the  general  problem.  For 
purposes  of  tractability  at  this  initial  stage,  scattering  from  elastic 
targets  is  considered  for  rotationally  symmetric  homogeneous  targets  in 
free  water.  The  problem  of  scattering  in  the  presence  of  a  bottom  is  solved 
for  sound-hard  and  -soft  targets  on  a  hard  or  soft  bottom. 

TV 

1.1  HISTORICAL  BACKGROUND  -  COMPARISON  OF  METHODS 

The  calculation  of  acoustic  scattering  from  an  elastic  solid 
submerged  in  a  fluid  needs  to  take  into  account  the  compressional  and  shear 
waves  inside  the  solid  as  well  as  the  compressional  wave  in  the  medium. 

The  problem  was  first  solved,  analytically,  in  the  frequency  domain,  for 
the  special  case  of  the  sphere  by  Faran  [1].  His  method  consists  of  writing 
the  general  eigenfunction  expansions  of  the  three  waves  in  the  frequency 
domain  and  solving  for  the  coefficients  by  applying  the  boundary  conditions. 
Hickling  [2]  published  results  of  calculations  based  on  Faran 1 s  solution  for 
spheres  of  several  different  materials.  An  eigenfunction  solution  is  re¬ 
stricted  for  use  with  very  special  geometries,  such  as  the  sphere  and  infi¬ 
nite  cylinder  at  normal  incidence  considered  by  Faran.  The  result  of  this 
method  is  a  complete  frequency  response.  Prior  to  Faran' s  work,  the  problem 
was  solved  in  a  similar  fashion  by  Anderson  [12]  for  fluid  targets.  Fluid 
targets  lack  a  shear  wave  and  can  be  considered  a  special  case  of  elastic 
target.  As  part  of  the  current  effort,  calculations  were  made  based  on  the 
above  solutions  for  an  elastic  sphere  in  order  to  check  out  the  solution  by 
the  more  general  formulation  to  be  described  below. 
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Another  approach  is  to  write  the  solutions  to  the  wave  equations 
as  integrals  over  the  bounding  surface  using  Green's  function.  This  method 
has  been  used  extensively  for  electromagnetic  scattering  problems  and  to 
some  extent  for  acoustic  problems.  For  scattering  from  a  perfect  conductor 
in  the  electromagnetic  case,  or  from  a  perfectly  sound-hard  or  sound-soft 
target  in  the  acoustic  case,  this  method  results  in  one  integral  equation. 
The  equation  is  solved  first  for  the  surface  field  values,  and  is  then 
used  with  surface  values  to  compute  the  scattered  field.  The  solution  of 
the  integral  equation  on  the  surface  is  usually  performed  in  the  frequency 
domain;  either  by  using  eigenfunction  expansions  [3]  (which  has  the  same 
geometry  restrictions  as  noted  above) ,  or  by  integrating  numerically  over 
the  surface  [4].  Numerical  integration  frees  the  problem  of  geometry 
restrictions.  The  frequency  domain  solution  requires  the  inversion  of  an 
NxN  matrix,  where  N  is  the  number  of  sample  points  on  the  surface,  for 
each  value  of  frequency. 

The  approach  taken  at  the  Sperry  Research  Center  is  to  solve 
the  integral  equation  numerically  directly  in  the  time  domain  [5] ,  by 
stepping  the  solution  in  time;  a  technique  named  the  space- time  integral 
equation  (STIE)  approach.  This  procedure  eliminates  the  need  for  a  matrix 
inversion.  The  time  domain  response  can  then  be  convolved  inversely  with 
the  excitation  to  obtain  the  frequency  response  of  the  target.  This  resul¬ 
tant  broadband  response  is  valid  from  zero  frequency  to  an  upper  limit 
determined  by  the  frequency  content  of  the  excitation  pulse.  The  latter 
is  limited  primarily  by  computation  time.  Practical  values  of  the  upper 
frequency  limit  are  ka = 6  to  12.  As  this  range  of  frequencies  includes 
the  resonance  region,  it  is  particularly  capable  of  relating  target  response 
to  size,  shape,  and  material.  The  relative  merits  of  frequency  and  time 
domain  approaches  are  discussed  in  [6] . 

For  the  elastic  problem,  three  integral  equations  are  required: 
two  scalar  equations  for  the  compressional  waves  in  the  target  and  medium, 
and  a  vector  equation  for  the  shear  wave  in  the  target.  These  equations 
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must  be  solved  simultaneously,  subject  to  the  boundary  conditions.  Depend¬ 
ing  on  the  variables  used,  these  can  lead  to  rather  complex  integral  equa¬ 
tions.  Hence,  the  first  problem  is  the  appropriate  choice  of  field  varia¬ 
bles.  In  this  report,  the  curl  and  divergence  of  the  displacement  are 
chosen,  since  these  yield  relatively  straightforward  integral  equations 
and  boundary  relations.  Fluid,  sound-hard,  and  sound-soft  targets  can  be 
treated  as  degenerate  cases  of  the  general  elastic  target  in  this  formula¬ 
tion. 


Other  approaches  of  interest  are  the  physical  optics  approximation 
and  the  geometrical  theory  of  diffraction  (GTD) .  The  physical  optics  ap¬ 
proximation  is  based  on  the  observation  that  in  the  limit  of  very  high 
frequency  (for  a  sound-hard  or  sound-soft  target) ,  the  ramp  response  of  a 
target  is  proportional  to  the  area  function  (i.e. ,  the  area  illuminated 
by  the  incident  ramp  as  a  function  of  time)  [7].  This  approximation  is 
used  here  for  predicting  the  asymptotic  behavior  of  the  response  as  the 
frequency  increases.  The  GTD  proceeds  by  calculating  the  response  and 
relative  phase  of  each  discontinuity  (scattering  center)  of  the  target  and 
then  combining  the  separate  responses  at  the  appropriate  relative  phases. 
This  theory  applies  only  at  high  frequencies  since  the  interaction  of  the 
surface  with  itself  (creep  waves)  is  neglected;  however  it  predicts  readily 
the  effect  of  aspect  change  on  the  response. 

1.2  VERIFICATION 

The  general  STIE  solution  has  been  checked  for  the  elastic  sphere 
by  comparison  with  the  eigenfunction  solution.  In  addition,  the  calculated 
responses  were  compared  with  measurements  taken  at  NCSL  on  a  variety  of 
targets.  The  results  will  be  discussed  in  Sections  4  and  5. 
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SECTION  2 

THEORETICAL  FOUNDATION 

2.1  PHYSICS  OF  ELASTIC  SOLIDS 

The  derivation  of  the  acoustic  wave  equations  in  fluids  and  elas¬ 
tic  solids  will  be  presented  in  outline  in  order  to  make  clear  the  meaning 
of  certain  variables  to  be  used  in  the  development  of  the  STIE  solution  in 
the  next  section.  The  reader  is  referred  to  a  modern  text,  such  as  [8]  for 
more  detail. 

2.1.1  Generalized  Hooke's  Law 

Consider  an  element  of  matter,  infinitesimal  but  continuous  (i.e., 
not  a  particle) ,  at  position  r  in  a  cartesian  coordinate  system.  Its  dis¬ 
placement  from  equilibrium  position  is  given  by 


u  =  dr  , 


and  the  total  deformation  by  du.  To  obtain  the  true  deformation,  or  strain, 
the  change  in  displacement  due  to  a  pure  rotation  about  r  is  subtracted 
from  du,  namely  ~(V x  u)  xu  [8].  Forming  the  total  vector  differential  and 
subtracting  the  rotation  term,  one  obtains  the  strain  tensor  D: 


D  D  D  1 

xx  xy  xz 

i 


9u 

.  /3u  3u  \ 

.  /3u 

3u  V  1 

_ X 

I  +  _x 

A  _£  . 

3x 

2  \  3  x  3y  / 

3u 

2  \  3x 

9z  /  I 

1 

I 

3 

• 

3u 


3z 


(1) 
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so  that  the  true  deformation  can  be  written  as  the  matrix  product 


du 

true 


[x,  y ,  z]  [D] 


dx  ~j 
dy  j 
.dz  J 


where  x,  y,  z  are  the  cartesian  unit  vectors,  and  ds  =  x  dx  +  y  dy  +  z  dz  is 
the  differential  change  of  position  within  the  infinitesimal  element.  The 
interpretation  of  the  diagonal  terms  of  D  is  that  Su^/Sx  is  the  linear 
(or  "dilational")  stress  component  in  the  x  direction,  i.e. ,  the  fractional 
increase  in  length.  The  off-diagonal  terms,  for  instance 

/3u  3u 

D  =  ihr*  +  ^ 

xy  2  \  dx  dy 


represents  one  half  of  the  sum  of  the  two  shearing  angles  about  the  z-axis 
(see  Fig.  2b) . 

The  stress  tensor  S  is  written 

S 

xz  i 

(2) 

S 

zz 

Here,  S  is  the  normal  stress  on  face  z  in  the  z  direction  (Fig.  2c) ; 
zz 

S  is  the  shear  stress  on  face  z  in  the  x  direction; 
zx 

S  is  the  shear  stress  on  face  z  in  the  y  direction, 
zy 

The  form  of  Hooke's  Law  is  that  Modulus  "times"  strain = stress.  For  the 
dilational  strain  the  physical  relation  is  (in  a  homogeneous,  isotropic  solid) 


S  = 


S  s 
xx  xy 


Y  D 


XX 


=  S 


XX 


C(S 


yy 


s  ) 

zz 


(3) 


i.e.,  the  linear  strain  in  the  x  direction  is  proportional  to  the  sum  of 


the  normal  stress  on  the  x  face  plus  a  contraction  proportional  to  normal 
stresses  applied  to  the  other  faces.  Y  is  Young's  modulus,  0  is  Poisson's 
ratio.  Similar  expressions  hold  for  the  other  diagonal  terms  of  D.  For 
the  shear  strains,  we  have  for  example 


2\i  D 


xy 


S 

xy 


(4) 


where 


p  is  the  shear  modulus, 

(the  factor  of  2  here,  is  due  to  the  factor  1/2  in  the  definition  of  D  ) . 

xy 

Equation  (3)  can  be  solved  for  S  ,  etc.,  so  that  (3)  and  (4)  can  be  written 
as 

S  =  A  (V  •  u)  I  +  2y  D  (5) 

This  is  the  generalized  Hooke's  Law  relating  stress  and  strain;  where 

/  y  y 

A  is  known  as  the  Lame  coefficient,  A  =  — — — — - — - -  , 

3(1-20)  3(1  +  0)  ' 

(V  •  u)  is  the  dilation  (or  total  linear  strain) . 

The  upper  diagonal  component,  for  example,  of  this  equation  is 

9u 

s*x '  +  211  ir  <5'> 


2.1.2  Wave  Equations  in  Elastic  Solids  and  in  Fluids 

The  x  component  of  the  total  force  on  element  (dx  dy  dz)  is  the  sum 
of  differences  between  opposing  faces  of  all  the  x  components  of  stress: 

/  9s  9s  9s  \ 

_ xx  _ yx  ,  _ zx  | 

\  9x  +  9y  +  9z  /dx  dy  dz  , 

with  similar  equations  for  the  y  and  z  components.  Using  equations  (5)  we 
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can  thus  write  in  vector  form 


Force/unit  volume  =  (A  +  y)  V(V  •  u)  +  y  u 


It  can  be  shown  [8]  that  for  small  ci  , 


where 


Force/unit  volume  =  p  u 


p  is  the  density  of  the  material. 


Equating  (6)  and  (7)  and  using  the  identity 


V  u  =  V(V  •  u)  -  VxVxu 


we  obtain  the  desired  equation  of  motion  in  a  solid: 


HH*) 


V  (V  •  u)  -  —  V  x  V  x  u 
P 


We  now  substitute 


2  A  +  2y  , 

=  — - -  (compressional) 


2  y 

c2  =  -  (shear) 


These  are  identified  as  the  compressional  and  shear  wave  velocities,  upon 
taking  the  divergence  and  curl  of  equation  (9) .  The  results  are  the  wave 
equations 


rl  9 

r~. 7  (V  •  u)  =  V^(V  •  u) 
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2  a  4.2 
c23t 


(V  x  u) 


=  V2  (V  x  u) 


(13) 


These  are  the  equations  to  be  solved  by  the  STIE  technique  of  Section  2.2 
They  state  that  the  divergence  and  curl  of  the  displacement  each  satisfy 
wave  equations.  In  the  fluid  medium  surrounding  the  solid,  there  is  no 
shear  stress,  so  that  we  have  only  equation  (12)  or  the  first  term  of  equ 
tion  (9) .  Using  3  as  a  subscript  denoting  the  medium: 


U3  =  P3  V(V‘G3)' 


or 


where 


-27-2  =  v2<?-u3)  . 

c  Jt 
3 


(14) 


(15) 


2.1.3  Elastic  Constants 

An  isotropic  solid  is  completely  characterized  by  two  elastic 
constants  and  its  density.  A  fluid  requires  one  constant  in  addition  to 
the  density.  For  convenience,  the  interrelationships  between  constants 
are  collected  here. 


Y  =  Young  1 s  modulus 
B  =  Bulk  modulus 
U  =  Shear  modulus 
0  =  Poisson's  ratio 
X  =  Lame  coefficient 


The  moduli  are  usually  given  in  units  of  10  psi.  A  convenient  conversion 

3  2  2 

factor  is  1  psi  =  6.8948  (g/cm  ) (m/s)  .  often  the  units  given  are  kg/mm  , 

where  kg  means  kg-force.  Such  values  are  to  be  multiplied  by  9.8  m/s2  to 

,  2  3  2 

obtain  the  correct  units,  i.e.  1  kg-force/mm  =  9800.  (g/cm  ) (m/s)  . 

Poisson's  ratio  is  dimensionless.  The  interrelationships  are 


Y  =  3B(1  -  20) , 

Y  =  2y  (1  +  0)  , 

X  =  B-2^  , 


2  X  +  2u  . 

c^  =  — - - &  (compressional  velocity)  , 


2  u 

c  =  -r  (shear  velocity) 
P1 

c2/c2  -  2(1-0) 
r  2  ~  (1  -  20)  ' 


0  = 


iff L 


For  a  fluid  one  can  consider  y  =  0  and  a  =  1/2.  Using  the  subscript  3  for 
the  fluid  medium: 


X 

c 


3 

2 

3 


=  B. 


B3 

=  — —  (compressional  velocity  in  fluid) . 
P3 


In  what  is  to  follow,  the  speeds  of  sound  and  densities,  p,  are  used  as 
the  elastic  constants. 


2.1.4  Boundary  Conditions 

The  boundary  conditions  at  the  interface  between  medium  and 
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and  target  are 


(1)  continuity  of  normal  displacement, 

(2)  continuity  of  normal  stress, 

(3)  continuity  of  tangent  shear  stress. 

These  conditions  are  very  general.  The  first  states  that  the  interface 
does  not  separate,  the  others  that  the  interface  does  not  contain  any 
stresses  of  its  own.  (The  latter  conditions  would  not  hold,  for  example, 
if  surface  tension  were  important,  or  if  the  boundary  were  a  thin  shell. 
In  these  cases,  the  stresses  in  the  boundary  surface  would  equal  the  dis¬ 
continuity  in  the  stresses  between  medium  and  target.) 

Equating  the  values  of  these  quantities  on  inside  and  outside: 


1.  u  =  u  , 
n3  n 


2.  P3c3  (V  •  u3)  =  P3c3  (V  •  u) 


2P1C2  (VG)  ' 


!.  o  =  y(- 


3n  3tanj 


Condition  (2)  is  obtained  from  the  normal  diagonal  component  of  Hooke's 
Law  (equation  5) ,  using  the  relation 


IT-  ,v‘5)  •  (Vs ' “> • 


The  operator  Vg  contains  only  surface  derivatives.  Condition  (3)  is  ob¬ 
tained  from  the  off-diagonal  terms  in  the  same  row  of  equation  (5) . 

Additional  boundary  relations  which  can  be  obtained  from  the 
above  and  the  equations  of  motions  are 


ii' »  n  wimwmppwwpw 


10 


-}■ 


SPERRY 

RESEARCH  CENTER 


cl  (V  •  u,)  =  cj  (V  •  u)  +  2  cjv?  u 
’  1  3n  2  S  n 


'3  3n 


5 .  (n  •  V  x  u)  =0 


where 


n  is  the  unit  normal  out  of  the  surface. 

Condition  (4)  is  obtained  by  equating  the  normal  components  of  the  equations 
of  motion  (9)  and  (14)  (justified  by  condition  (1)),  and  by  using  the  rela¬ 
tion 


n  •  V  x  (V  x  u)  =  n  •  V  x  (V  x  u)  =  -V  •  n  x  (V  x  u)  =  -2  V^u  .  (16) 

S  S  S  n 

The  last  step  of  the  above  is  from  the  fact  that  boundary  condition  (3) 

implies  that  the  surface  components  of  the  curl  can  be  written  in  terms  of 

surface  derivatives  of  u  : 

n 

n  x  (V  x  u)  =  2  V  u  .  (17) 

S  n 

Condition  (5)  can  be  derived  through  use  of  condition  (3). 

It  is  to  be  noted  that  the  boundary  conditions  also  hold  under 

differentiation  by  t  and  directional  differentiation  in  the  surface  (V  ) . 

s 

Also  note  that  by  equation  (14)  and  B.C.  (1): 


(V  •  uj 


(18) 
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For  a  fluid  target  these  boundary  conditions  are  simplified  by 
the  facts  that  =  0  and  that  condition  (3)  does  not  apply  (or  rather, 
is  satisfied  trivially  by  the  fact  that  y  =  0)  .  We  thus  have,  for  a  fluid 
target  in  a  fluid  medium: 


1.  u  =  u 
n3  n 


2.  P3C3  *  U3J  =  pici  *  “> 


4-  C3  ^  (V  *  V  =  (V  ’  S) 


In  addition,  (Vxu)  is  identically  zero  in  a  fluid. 


As  a  further  degeneracy,  consider  sound-hard  and  soft  targets, 
for  which  no  acoustic  energy  penetrates  the  target,  i.e.,  V  •  u  is  also 
identically  zero.  A  sound-hard  target  can  be  considered  as  the  extreme 
case  where  p^  -►  hence  we  obtain  the  boundary  condition 


u 

n 


2  9 

°3  IS  <V 


»3>  ■ 


and  only  (V  •  u^)  is  free  to  vary.  A  sound-soft  target  can  be  considered 
as  the  other  extreme  where  p  0,  so  that 

V  •  u3  =  0 

and  only  u  is  free  to  vary. 
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2.1.5  Relations  Between  Displacement  and  Pressure 

In  an  ideal  fluid  there  are  no  shear  stresses  and  Hooke's  Law 

reduces  to 

P6  =  -B (V  •  u3)  , 

or  (19) 

Pe  -  -Cj  P3  (V  •  u3> , 
e 

where  P  is  the  excess  pressure.  The  minus  sign  is  introduced  since  (V  •  u) 
is  the  change  in  volume/unit  volume,  which  is  negative  for  an  increase  in 
pressure. 


At  the  surface 

3pe  2  3  „  _ 

17"  ‘°3  P3  to  <V'"3>  *  -p3“n-  <20) 

Hence  in  the  fluid,  the  displacement  can  be  interpreted  directly  in  terms 
of  excess  pressure. 

2.2  INTEGRAL  EQUATION  SOLUTION 

2.2.1  STIE  For  Scattered  and  Interior  Waves 

Using  the  Green’s  function  G(r|r',  t|x)  =  —  6 (t  -  R/c  -  t) , 

_  R  3 

where  R  =  r  -  r  ,  a  solution  to  the  wave  equation  can  be  developed,  known 
as  the  Kirchoff  integral  representation: 


V  •  u3(r,t) 


fin'-hL 


3<V5'>dS' 


(21) 
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where  the  integrals  are  over  the  bounding  surface, 

n^  is  the  unit  normal  out  of  the  surface, 

R  is  the  magnitude  of  R  =  r-r'  , 

A  _ 

R  is  the  unit  direction  of  R  , 

t  =  _L  +  _L_  _L 

3  ~  R2  C3R  8T  ' 

(V  •  U3)  =  V  •  u3 (r',T)  ,  with  T  =  t  -  R/c3  . 

inc 

(V  •  u3>  is  due  to  sources  external  to  the  surface 
(the  incident  pressure) . 

The  above  is  for  the  exterior  wave.  Maintaining  the  convention  that  n 
points  out  of  the  surface,  for  the  interior  waves,  the  sign  of  the  inte¬ 
grals  is  changed  and  there  is  no  incident  term.  Hence: 


V  •  u(r,t) 


•-  f  1 

41T  J  R 


—  (V  •  u^)  dS ^ 


3n 


-  —  f 

4tt  J 


(n  •  R)  Z,  (V  •  u')dS' 


(22) 


and 


V  x  u(r,t) 


4iT  J  R 


1  ^  <  r~i  ™  ^ > 

9^(Vxu)dS  4TT 


5  /  ir 


R)  L2(Vxu')dS'  (23) 


Here,  L L^,  (V*  u^) ,  etc.  are  defined  analogously  to  the  quantities  in 
equation  (21).  Note  that  the  retardation  is  t-T  =  R/c^  in  equation  (22) 


and  t-T  =  R/c2  in  equation  (23). 


2.2.2  ST IE  on  Boundary 


To  solve  the  scattering  problem,  these  equations  are  first  solved 
on  the  boundary,  after  which  the  scattered  field  may  be  obtained  from 
equation  (21).  To  obtain  the  solution  on  the  surface,  the  observation 
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point,  r,  is  moved  to  the  surface  and  a  limiting  process  applied.  The 
result  is  the  following  set  of  integral  equations  for  r  on  the  surface: 


V'Vr't)  =  2(v.u3(?,t))lnc  -  Ay*  I_i,  (V  •  u3)ds" 

+  y  (n”"  •  R)  L3(V*u3)dS^  (2- 

v  •  G<;,t,  .  i  y  I  ^  .  5-lds-  .  J-  y  (8- .  6,  £i(v  •  3')ds'  (2i 

V  x  u(r,t)  =  A y  i  (V  x  u")dS"  -  A  y  (ft"  .  R)  x  u')ds'  (2( 


To  solve  the  above,  the  integrals  are  replaced  by  summations 
over  finite  patches,  over  which  R  and  the  field  variable  are  assumed  con¬ 
stant.  The  "self-patch"  contribution  (which  includes  the  point  r)  is  in¬ 
tegrated  analytically.  For  the  first  term  of  the  integral  over  the  self¬ 
patch  at  r^,  we  define 


For  the  second  term,  we  represent  V  •  u3  by  a  Taylor  expansion  about  r: 
V  •  u3  =  v  .  53  +  p(A  (7  •  u3>  cos  e  +  A  (V  .  u3)  sin  e  -  -^A_  (V  .  u3)j 

and  n'  by 


n'  =  n  -  p (K  x  cos  0  +  K  y  sin  0)  , 


where 

Kl'  K2  are  Pri-nciPal  curvatures, 

x,  y  is  the  surface  (curvilinear)  coordinate  system,  dS^  =  pdpd8. 

The  second  integral  over  the  self-patch  then  becomes  e  (V  •  u  )  • 

k  3 

where 


(28) 


We  then  can  write 


(1  -  S)  (V  •  53,k  ♦  Y*  £  <v  •  53,k  -  -I3k  *  2,7  •  u/™ 


(1  +  0  (V  •  u).  -  Yu  f-  (V  •  u)  =  I., 
k  k  k  dn  k  lk 


(29) 


,l*ek)(7xG)k-Yk-,7x=)k-I2k 


where- 


(V 


u_>  An.  •  R.J 
33  3  3k  I 


with  similar  expressions  for  1^  and  I2k  (the  "non-self  integrals") . 

The  discrete  formulations  of  the  non-self  integrals  contain  ten 
quantities:  V  •  u^,  V  •  u,  V  xu  and  their  normal  derivatives  at  retarded 
times.  One  quantity  is  identically  zero  (B.C.  5) ,  so  that  9  quantities 
need  to  be  stored  at  all  space  and  time  sample  points  to  evaluate  the  in¬ 
tegrals.  If  the  time  step  is  chosen  sufficiently  small  so  that  the  retar¬ 
dation  is  always  greater  than  a  time  step,  the  solution  can  proceed  by  a 
straightforward  process  of  marching  in  time. 
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2.2.3  Simultaneous  Solution  of  STIE  on  Boundary 


At  each  space-time  point,  equations  (29)  need  to  be  solved  simul¬ 
taneously  for  the  9  quantities  of  interest.  This  is  done  by  application 
of  the  boundary  conditions.  Using  B.c.  2  and  4,  the  first  two  equations 
can  be  rewritten  (suppressing  the  subscript  k) : 


(1-e)  (V  •  u  ,)  +  Y  ~~  (V  •  uj  =  2 (V  •  u)inc  -  I, 
3  on  3  3  3 


0(1  +  e)  (V  •  u3)  -  Y  (V  •  u3)  =  i1  -  a  |(l  +  e)  (Vg  •  u)  +  yV^  u^) 


(30) 


Here  we  have  used  the  relations 


(V  •  u)  =  B(V  •  u  )  +  a(V  •  u) 
3  s 


3n 


(V  •  u) 


-f£  (V 

C1 


V 


-  a  V" 


u  , 
n 


(31) 


where 
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This  pair  of  equations  can  be  solved  for  (V  •  u^)  and  its  normal  derivative, 
if  the  displacement  vector  u  is  known  at  time  t  over  the  entire  surface. 
Since  from  (18) , 

a2 

3  Un  Ii(n 
a.  2  2  9n  (V  U3)  ' 

3t  c3 

u  can  be  found  from  earlier  values  by  the  relation 
n  2 

a2 

2  3  u 

u  (t)  =  2  u  (t-At)  -  u  (t  -  2  At)  +  (At)  - -  (t-At).  (32) 

n  n  n  9t2 

The  surface  components  of  u  can  be  found  by  numerical  solution  of  the 
interior  wave  equation  (9).  For  example,  if  ^  and  a are  orthogonal 
directions  in  the  surface  and  u^  is  a^  component  of  u,  etc: 


3^ 

=Cl^(V*5)+C2  [h  (VxII)a2  -  9^  {Vx",n]*  (33) 

Since  9/9n  (Vxu)^  =  0,  the  normal  component  of  the  last  of  equations  (29) 
yields  directly 


(VxS,n  =  lT7Tr  (vfi)  • 


Surface  differentiation  then  yields  the  third  term  in  equation  (33).  Now 
the  surface  components  of  (V  x u)  can  be  obtained,  by  virtue  of  B.C.  3, 
from  equation  (17).  Hence  the  surface  components  of  the  last  of  equations 
(29)  yields 


9  t n  (1  +  e)  1  .=  .  ~  . 

5n  (,*U)a2  ■  y  IVxU,a2  *  Y  H2  *2>  ' 
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with  a  similar  expression  for  the  a^  component.  Hence  and  u^,  in  the 
surface,  can  be  found  by  equation  (33)  from  earlier  values  in  time  by 
application  of  a  relation  similar  to  equation  (32) . 

The  solution  is  now  complete.  It  can  be  seen  to  require  three 
integrations  over  the  retarded  field  values  on  the  boundary,  one  integra¬ 
tion  for  each  kind  of  wave  (exterior,  interior  compressional ,  interior 
shear) ,  where  each  wave  is  associated  with  its  own  retardation.  It  also 
appears  that  9  quantities  need  to  be  stored  at  all  time  and  space  points. 
Although  these  quantities  are  related  by  the  boundary  conditions,  it  is 
computationally  advantageous  to  store  them  all,  rather  then  recomputing 
them  when  required  to  evaluate  the  STIE. 


2.2.4  Incident  Field 

\ 

The  incident  pressure  is  taken  as  a  plane  wave  with  time  depen¬ 
dence 


a  -(at)' 
(V  •  u)lnC  =  —  e  n 

/if 


(36) 


the  "smoothed  impulse".  The  calculation  is  started  at  a  time  such  that 
the  incident  field  is  still  very  small  at  the  closest  point:  (a^t)-^. 

In  the  frequency  domain,  the  smoothed  impulse  has  the  same 

shape,  with  appreciable  amplitude  out  to  k  !U  tt  a  .  The  frequency  response 

of  the  target  could  thus  be  obtained,  by  convolution,  from  0  to  k  IT  a  . 

n 

Since  larger  a  requires  smaller  time  and  space  steps  for  solution  accuracy, 
n 

practical  considerations  (computer  storage  and  running  time) ,  limit  a 

n 

to  1  or  2. 

For  the  computations  of  a  target  on  the  bottom  (see  Section  2.4) 
a  nearby  source  was  taken,  which  includes  a  1/R  dependence.  The  amplitude 
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of  this  spherical  wave  was  normalized  such  that  if  the  source  were  located 
at  00  on  a  hard  bottom,  the  result  would  be  the  incident  field  of  equation 
(36) .  For  nearby  source: 


.me 


a  R 
n  s 


~(a  t  )' 
n  1 


(at,)' 

n  2 


/if  ! 


(37) 


where 


+  is  for  the  hard  bottom, 
-  is  for  the  soft  bottom. 


Rg  =  distance  from  source  to  target  center. 


R^  =  distance  from  source  to  target  point, 

R2  =  distance  from  image  source  to  target  point. 


tl  =  t  +  Rs-Rr 


t  =  t  +  R  -  R  . 
2  S  2 


2.2.5  Scattered  Field 


The  scattered  pressure,  P  ,  or  equivalently  (V  •  u^) ,  is  given  by 


equation  (21)  of  Section  2.2.1;  repeated  here 


V  •  u3(r,t) 


-  f~ 

4tt  J  r 


it <?  •  =? as'  + 


i 


•  R)£3<V  •  u3)ds"  (38) 


which  is  valid  in  the  near  and  far  field.  That  is,  of  all  the  quantities 

—  2  2 

computed  on  the  surface,  only  two:  (V  •  u,)  and  9/8n  (V  •  u,)  =  1/c  9  u  /c 

3  3  3  n 

need  to  be  retained  for  the  scattered  pressure  calculation. 
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In  the  limit  as  r  becomes  large  the  above  becomes 


32 

r(v  •  u3(I.t))  =  i  /  -  i  -£?  *  <«'  ’  f  T  £  ^ 


(39) 


Note  that  the  (bistatic)  response  can  be  calculated  for  any  observation 
direction  without  recalculating  the  surface  quantities;  the  latter  are 
associated  with  a  particular  direction  of  incidence.  The  near  field  eval¬ 
uation  of  (38)  was  done  for  the  case  of  targets  on  the  bottom. 

2.3  REDUCTION  OF  ELASTIC  SCATTERING  PROBLEM  TO  SIMPLER  CASES 

It  is  instructive  to  see  how  these  equations  reduce  to  simpler 
forms  for  simpler  problems. 

2.3.1  Fluid  Target 

For  example,  consider  the  fluid-fluid  problem,  for  which  the 

target  does  not  support  shear  waves.  For  this  case,  the  constants  y  and 

c2  are  zero,  and  there  is  no  equation  (23).  Boundary  condition  (3)  does 

not  apply  and  boundary  conditions  (2)  and  (4)  are  simplified  by  the  loss 

_  2 

of  exactly  those  terms  (V„  *  u  and  V  u  )  which  required  the  numerical  solu- 

S  S  n 

tion  of  the  wave  equations  in  the  elastic  problem.  Equation  (29)  can  thus 
be  solved  directly,  using  boundary  relations  (2)  and  (4).  The  retarded 
integration  requires  storage  of  only  2  quantities  in  all  space  time.  Also, 
using  relation  (19) ,  the  equations  can  be  written  in  terms  of  the  excess 
pressure  P  and  its  normal  derivative  9P/9n.  Equation  (24)  for  pressure 
on  the  surface,  then  becomes 


P(r,t)  =  2PinC  (r 


•°  -  w/r 


__9_ 

9n' 


P  dS 


—  f 

2-nJ 


L3  p"  (n'  •  R)ds'  . 


(40) 


Equation  (25)  can  be  similarly  reformulated. 
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2.3.2  Hard  and  Soft  Targets 

If  no  acoustic  energy  penetrates  the  targets,  equation  (25)  also 
disappears  and  we  are  left  simply  with  (21)  or  (40) .  The  boundary  condi¬ 
tions  for  this  case  are  either  that  the  target  is  perfectly  rigid  (sound- 
hard)  for  which  9p/9n=0,  or  that  the  target  is  a  pressure-release  surface 
(sound-soft)  at  which  p  =  0.  For  both  cases,  one  of  the  integrals  in  (40) 
disappears.  The  solution  thus  requires  only  one  space-time  integration 
over  the  retarded  boundary  value  and  requires  the  storage  of  only  one 
quantity  in  space- time.  This  is  the  problem  which  has  been  solved  at  the 
Sperry  Research  Center  for  a  variety  of  geometrical  configurations  [9,10,11] 


2.4  TARGET  IN  THE  PRESENCE  OF  A  BOTTOM 


We  consider  here  the  case  of  perfectly  hard  or  soft  targets;  that 
is,  equation  (40)  applies  with  the  condition  that  either  9p/9n =  0  (sound- 
soft)  or  p =  0  (sound  hard).  After  solving  analytically  for  the  "self-patch" 
this  equation  is  rewritten  as 


P (r , t) 


1 

(1  -  e) 


2  plnC  (r,t) 


(41) 


for  the  sound-hard  target;  or  as 


!e 

9n 


(r,t) 


1 

Y 


~  inc 

2  p 


(r,t) 


_1_  f  I  9  - 

2TT  J  R  9n^  p 

NS 


(42) 


for  the  sound-soft  target,  where  the  integrals  are  over  the  non-self  part 
of  the  surface  and  y,  e,  L are  defined  in  Section  2.2.2. 

If  a  sound-hard  or  soft  plane  near  the  target  is  added  to  the 
scattering  problem,  rather  than  solving  (41)  or  (42)  on  the  plane,  a  method 
of  images  is  applied. 
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2.4.1  Method  of  Images 

We  assume  again  that  the  plane  bottom  is  either  perfectly  hard  or 
perfectly  soft.  That  is,  we  have  one  of  equations  (41)  or  (42)  on  the  bot¬ 
tom.  Also  assuming  that  the  plane  is  00 ,  we  see  in  Fig.  3  that  the  boundary 
conditions  on  the  plane  are  satisfied  if  we  remove  the  bottom  and  add  sym¬ 
metrically  placed  targets  and  sources  below  the  bottom  plane.  This  latter 
configuration  has  the  same  pressure  field  in  the  upper  half-space  as  does 
the  original  configuration.  In  the  case  of  a  hard  bottom,  the  reflected 
source  has  the  same  sign  as  the  original  source;  and  in  the  case  of  a  soft 
bottom,  the  reflected  source  must  be  of  opposite  sign  to  satisfy  the  bound¬ 
ary  condition  on  the  plane. 

We  note  that  by  symmetry  the  value  of  the  pressure  (or  its  normal 
derivative)  at  a  point  on  the  image  target  is  the  same  as  that  at  the  sym¬ 
metric  point  on  the  original  target  (with  a  +  or  -  sign,  respectively,  for 
a  hard  or  soft  bottom) .  In  the  numerical  calculation  we  therefore  do  not 
calculate  the  image  pressure  explicitly,  although  its  contribution  to  the 
pressure  at  a  target  point  is,  of  course,  included  in  the  summation  of 

equal  ions  (41)  and  (42) .  Also,  in  equations  (41)  and  (42) ,  we  understand 
inc 

p  to  mean  the  total  incident  pressure  wave  from  the  source  and  its  image 
(see  Section  2.2.4).  We  can  then  apply  these  equations  without  further 
change. 

Note  that  in  the  case  of  a  sound-soft  bottom,  the  effect  of  the 
source  and  its  image  tend  to  cancel  at  a  point  on  the  target.  The  pressure 
at  the  target  is  due  to  the  difference  in  path  length  from  the  two  sources  - 
the  closer  the  source  is  to  the  bottom,  the  smaller  is  the  target  surface 
pressure  until  the  source  hits  the  bottom,  which  is  equivalent  to  no  source 
at  all  and  the  pressure  is  zero  everywhere. 

In  the  case  of  a  sound-hard  bottom,  an  interesting  case  is  that 
for  which  a  target  (symmetrically  shaped)  is  half-penetrating  the  bottom 
with  the  source  located  on  the  bottom  (see  Fig.  3c) .  This  case  is  equivalent 
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FIG.  3  Method  of  images. 


to  the  target  in  free  space  with  double  the  strength  of  the  source.  This 
configuration  (with  the  source  moved  far  away)  was  used  to  check  the  new 
program  ECHOB  by  comparison  with  the  existing  program  ECHO  for  free  targets. 

2.4.2  Self-Term  Integration  for  Patches  Close  to  Bottom 

As  has  been  discussed,  the  non-self  integrals  in  the  STIE  are 
evaluated  numerically  by  dividing  the  target  into  a  number  of  patches  and 
assuming  constant  values  of  R  and  p  or  3p/9n  over  a  patch.  This  becomes 
inaccurate  when  R  is  small,  which  occurs  in  particular  when  the  integration 
point,  r' ,  is  the  image  of  the  observation  point,  r,  and  r  is  close  to  the 
bottom.  This  happens,  for  example,  when  a  sphere  or  cylinder  lies  tangen¬ 
tially  to  the  bottom.  This  problem  is  handled  by  integrating  such  patches 
analytically.  The  result  of  the  image  patch  integration  is  included  in 
the  expression  for  y  or  £  and  the  integrals  (41)  and  (42)  are  understood 
to  omit  this  patch.  The  result  depends  on  whether  the  bottom  is  hard  or 
soft. 


For  a  hard  target  on  a  hard  bottom: 


(43a) 


for  a  hard  target  on  a  soft  bottom: 
es  =  e(l-/l+ (2z/Y)2  +  (2z/y))  +  (nz~  z^~1~  i1_~ 


(2z/Y) 


\ 


A  +  (2z/y)2  / 


(43b) 


For  a  soft  target  on  a  hard  bottom: 


•A1 


Yh  =  Y 1 1  +  /l+  (2z/y)  -  (2z/y)  I 


(44a) 


for  a  soft  target  on  a  soft  bottom: 


=  y(i  -  <£+ 


(2z/Y)  +  (2z/Y) 


(44b) 


where  z  is  the  distance  of  the  patch  center  from  the  bottom,  n z  is  the  com¬ 
ponent  of  the  target  normal  in  the  direction  of  negative  bottom  normal, 

K^,  are  the  principal  curvatures  of  the  patch.  These  expressions  are 
derived  in  reference  [10] . 
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SECTION  3 

NUMERICAL  SOLUTION 

3.1  REPRESENTATION  OF  THE  TARGET 

The  surface  of  the  target  is  represented  by  a  set  of  more  or  less 
uniformly  spaced  nodes.  Each  node  is  associated  with  an  area  ("patch")  and 
is  defined  by  position  r,  normal  n,  principal  curvatures  K^,  ,  and  direc¬ 

tions  of  principal  curvature  a^,  a^,  where  n'a^xa^.  For  a  target  rota- 
tionally  symmetric  about  the  z-axis,  a^  is  taken  as  lying  along  the  gen¬ 
erating  curve  (or  meridian)  and  a^  =  $  (Fig.  4).  These  quantities  are 
pregenerated  and  stored  for  each  node  by  the  geometry  section  of  the  compu¬ 
ter  program,  (with  the  vectors  represented  in  the  cartesian  coordinate 
system) . 


The  target  is  centered  at  the  origin  and  the  incident  field 
direction  is  taken  to  lie  in  the  z-y  plane.  A  larger  number  of  nodes 

increases  the  accuracy  of  solution,  but  is  limited  by  computation  time 
2 

(h*  N  )  and  computer  storage  requirements  (^  N) .  Node  spacings  as  large  as 
0.4  (or  ^  90  nodes  on  a  sphere)  have  been  found  to  be  adequate  for  uniform 
meshes.  For  long  thin  targets  (such  as  the  cylinder)  a  non-uniform  mesh 
can  be  used,  with  higher  node  concentration  at  regions  of  high  curvature 
(or  near  edges)  where  the  field  is  expected  to  change  most  rapidly. 

3.2  NUMERICAL  REPRESENTATION  OF  THE  STIE 

In  Section  2.2.2  the  STIE's  were  integrated  analytically  over 

the  self-patch.  For  discrete  time  and  space  points  this  is  written  (at 

time  t.  at  point  r  ): 

D  K 

d-Vf^'V)*  +  **  ^(V-=3ltj>)k  •  2(7- VV)inC  -  I3SSk  U5) 
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where 


3NS, 


2tt  SjiL-  3n  (V  *  u3(Tjki)K  "  (fik£*  V  L3^  *  U3 (T jki?  ' 

l*k  1  K*  ; 


(46) 


t.,„  =  t. 
jkl  j 


\JI/C3  ' 


L  =-L-+  1 


2  '  \ic3  9T 


Similar  forms  are  written  for  (V  •  u)  and  (V  x  u)  . 

If  the  time  step.  At,  is  taken  sufficiently  small  so  that 
>  At  for  all  (k£) ,  then  the  field  values  in  the  non-self  integrals 
are  strictly  retarded.  Hence  the  non-self  integrals  can  be  evaluated  at 
time  t_.  from  the  field  values  at  earlier  times.  It  is  thus  necessary  to 
store  the  nine  quantities  that  appear  in  the  integrals  for  each  node  and 
for  all  time  points  from  the  current  time  back  to  the  time  corresponding 
to  the  greatest  retardation.  For  efficient  use  of  computer  storage  this 
required  time  window  is  stepped  along  with  the  stepping  in  t^.  The  size 
of  the  time  window  is  given  by 


T  =  R  /c  . 

W  max  min 


(47) 


where  R  is  the  maximum  distance  between  any  two  nodes  and  c  .  is  the 
max  min 

minimum  of  c^,  c2»  c^.  The  evaluation  of  the  field  quantities  at  time  T 
is  done  by  interpolation  on  the  discreet  values. 

The  non-self  integral  for  the  curl  equation  is  a  vector:  Each 
contribution  to  the  sum  is  evaluated  in  the  local  coordinate  system  of 
the  integration  point,  £,  then  transformed  to  the  coordinate  system  of 
the  observation  point,  k,  to  permit  the  vector  addition  of  the  contribu¬ 
tions.  The  values  of  the  displacement  vector  u  at  all  nodes  is  also  stored 
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for  the  current  and  previous  time  point. 

The  solution  procedure  can  now  be  summarized:  at  time  t, 

(1)  Evaluate  the  three  non-self  integrals  (2  scalar  and  one 
vector  equation) ,  using  the  retarded  values  of  the  nine 
field  quantities,  (V  •  u^) ,  (V  •  u)  ,  (Vxu),  and  their  normal 
derivatives. 

(2)  Evaluate  (Vsun^  and  ^ s  *  by  numerical  differentiation 
in  the  surface,  using  the  current  value  of  u. 

(3)  Solve  the  first  two  equations  of  (29)  for  (V  •  u^)  and 
3/3n  (V  •  u  )  at  the  current  time  point  and  find  (V  •  u) 
and  3/3n  (V  •  u)  by  the  relations  (31) . 

(4)  Evaluate  the  normal  component  of  (V x u)  from  the  normal 
part  of  the  vector  STIE  (equation  34) . 

(5)  Evaluate  the  surface  components  of  (Vxu)  by  numerical  dif¬ 
ferentiation  of  u  in  the  surface,  using  equation  (17) . 

n 

(6)  Use  the  relation  (35)  to  find  the  surface  components  of 
3/3n  (Vxu).  (The  normal  component  is  identically  0.) 

(7)  Use  the  result  of  (3)  and  (4)  to  find  the  surface  deriva¬ 
tives  of  (V  •  u)  and  (Vxu)  . 

n 

(8)  From  the  relation  (33)  evaluate  the  surface  components  of 

2_  2 
3  u/3t  . 

(9)  Obtain  u  at  the  next  time  point  (t  +  At) ,  by  numerical 
extrapolation  (equation  32) . 

The  process  is  now  repeated  for  the  next  time  point. 

Strictly  speaking,  since  the  nine  field  quantities  are 'related 
by  the  boundary  conditions,  not  all  need  to  be  stored  to  evaluate  the 
non-self  integrals.  In  fact,  it  can  be  argued  that  the  problem  is  fully 
described  by  only  4  quantities  (and  their  derivatives) .  It  was  found 
to  be  more  efficient  and  accurate,  however,  to  retain  all  9  field  variables. 
An  alternate  approach  is  described  in  the  Appendix. 
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3.3  SURFACE  DERIVATIVES 


An  important  step  in  the  above  procedure  is  the  evaluation  of 
derivatives  in  the  surface.  It  should  be  noted  that  these  derivatives  do 
not  appear  in  the  degenerate  cases  of  fluid,  hard,  and  soft  targets,  so 
that  the  elastic  problem  represents  a  large  increase  in  complexity,  requiring 
the  evaluation  of  a  vector  integral  as  well  as  surface  differentiation. 


A  convenient  way  to  form  these  derivatives  is  to  consider  a  plane 
tangent  to  the  surface  at  the  node  k  and  then  projecting  the  adjacent  nodes 
and  the  vector  u  at  those  nodes  upon  the  plane.  A  second  order  Taylor  ex¬ 
pansion  of  u  is  written  about  node  k,  e.g.  for  u  : 


9u 


3u 


s2 
d  u 


U  (£>  =  u  (k)  +  T-5-  (Aa_)  +  3-^  (Aa  )  +  - f 

n  n  da,  1  3a„  2  ~  2 

1  2  3a, 


(AajT 


32u  (Aa_) ' 

_ n  _ 2 


3aT 


g2u 

+  3^2(Aai>  (Aa2> 


(48) 


As  there  are  generally  6  adjacent  nodes  in  addition  to  the  node  k  itself, 
the  above  can  be  solved  for  the  derivatives  by  a  least  square  fit.  This 
fitting  process  needs  to  be  performed  only  once  and  the  result  put  into 
the  form  of  a  6x7  matrix,  S.  The  derivatives  are  then  evaluated  when  re¬ 
quired  by 


tD]  =  S[un]  (49) 

where  [D]  is  the  column  of  derivative  values  at  node  k  and  [u  ]  is  the 

n 

column  of  values  of  u^  at  k  and  the  adjacent  nodes.  One  matrix  S  and  the 
indices  of  the  adjacent  nodes  are  stored  for  each  node. 

The  advantage  of  this  process  is  that  the  evaluation  of  the  STIE 
becomes  independent  of  the  geometrical  representation  of  the  target.  In 
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particular,  nodes  do  not  have  to  be  placed  along  coordinate  lines  as  re¬ 
quired  by  the  usual  finite  difference  techniques. 

For  problems  with  special  geometries,  such  as  a  sphere,  the 
derivatives  may  be  written  in  the  curvilinear  coordinate  system  and  eval¬ 
uated  using  their  finite  difference  representations,  rather  than  using 
the  plain  projection  method  above.  For  the  first  order  derivatives  the 
two  approaches  give  similar  results;  the  second  order  derivatives,  however, 
can  be  different.  As  discussed  later,  at  the  present  time  the  problem  of 
accurate  determination  of  second  surface  derivates  has  not  been  resolved 
and  appears  to  be  the  cause  of  instabilities  in  the  solution  procedure. 

3.4  SYMMETRY  RELATIONS 

To  increase  the  efficiency  of  the  solution,  symmetry  conditions 
are  used  whenever  possible.  Scattering  at  arbitrary  angles  of  incidence 
from  a  rotationally  symmetric  target  has  plane  symmetry  (Fig.  4) .  The 
symmetry  plane  is  taken  as  the  x=0  plane.  The  field  quantities  are  cal¬ 
culated  and  stored  only  for  nodes  at  x>0  and  their  values  at  x  <  0  inferred 
from  the  symmetry  conditions. 

The  values  of  the  scalar  quantities  (V  •  u) ,  etc.  at  (-x)  are  the 
same  as  at  (x) .  For  the  vector  quantities  the  following  table  summarizes 
the  symmetry  conditions,  where  -  indicates  a  sign  reversal  upon  reflection 
in  the  x  =  0  plane 


V  x  u  -  +  -  + 
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[Recall  (Fig.  4)  that  a^  is  identified  with  <f>,  a 2  is  in  the  direction  of 
increasing  z  and  n  is  the  outward  normal.] 

For  a  rotationally  symmetric  problem,  i.e.  axial  incidence  upon 
a  rotationally  symmetric  target,  the  symmetry  conditions  are  even  more 
restrictive.  We  then  have  that  the  scalar  quantities  are  constant  under 
rotation  about  the  axis  of  symmetry.  The  vector  quantities,  in  our  curvi¬ 
linear  coordinate  system,  rotate  along  with  the  target  rotation,  hence  the 
items  marked  in  the  table  on  the  preceding  page  are  required  to  be  0  for 
the  rotationally  symmetric  problem: 

Vector 


(V  x  u)  j 

3.5  INSTABILITY  PROBLEMS 

Most  of  the  instability  problems  that  arise  in  the  STIE  solution 
have  been  solved;  one  problem,  that  of  a  stable  formulation  of  surface 
derivatives,  remains. 

‘  3.5.1  Types  of  Instabilities 

Three  classes  of  instabilities  have  been  identified,  which  may 
be  named,  suggestive  of  their  origin:  time-,  space-,  and  resonance  insta¬ 
bilities. 

The  time- like  instability  is  affected  strongly  by  the  time  step 
of  the  solution,  setting  in  rapidly  when  the  time  step  is  too  large. 

The  onset  of  this  instability  is  made  later  by  taking  a  time  step  smaller 
than  the  time  required  for  a  wave  to  propagate  between  adjacent  nodes, 
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At  <  R  .  /c 

min  max 


(52) 


Reducing  the  time  step  further  does  not  affect  the  stability  appreciably. 

The  nature  of  this  instability  is  rapid  oscillations  which  build  up  in 
amplitude  very  rapidly  and  are  in  phase  on  the  surface.  The  calculated 
far  field  response  is  smooth  up  to  some  point  near  the  end  when  it  suddenly 
breaks  up  with  rapidly  increasing  oscillations.  This  problem  does  not 
occur  for  the  perfectly  hard  target,  probably  because  the  nearest  patches 
do  not  contribute  heavily  to  the  non-self  integral  (since  there  is  no 
9p/3n  term.  This  problem  is  completely  eliminated  by  including  some 
smoothing  in  the  time  interpolation  performed  in  the  evaluation  of  the  re¬ 
tarded  contributions_.to  the  non-self  integral.  The  time  interpolation  is 
performed  to  second  order  on  the  values  at  the  three  time  points  closest 
to  the  time  T ;  these  three  values  however  are  obtained  by  linear  smoothing 
of  the  nearest  five  time  points.  This  process  eliminates  the  instabilities 
as  they  arise  and  also  yields  continuous  first  derivatives  as  required. 
(Linear  interpolation  gives  rise  to  another  instability  due  to  discontinuity 
in  the  first  derivates;  higher  order  interpolation  does  not  provide  suffi¬ 
cient  smoothing . ) 

The  space-like  instability  appears  to  be  caused  by  the  discrete 
representation  of  the  target  geometry.  It  appears  as  an  inaccuracy  at 
some  critical  point  (the  specular  point  or  edges)  and  is  then  propagated 
with  slowly  increasing  amplitude  over  the  surface.  Oscillations  on  the 
surface  are  in  phase.)  The  period  of  oscillation  is  slower  and  corresponds 
to  the  target  size.  The  size  of  the  time  step  only  affects  the  phase 
but  not  the  magnitude  of  this  instability.  This  problem  can  be  eliminated 
by  modeling  the  target  with  a  finer  mesh  near  the  critical  points.  A  finer 
mesh,  by  relation  (52) ,  also  requires  a  smaller  time  step,  which  may  not 
be  practical  for  large  targets.  To  keep  the  time-step  large  (i.e.  only 
small  enough  to  satisfy  (52)  for  most  nodes)  an  iterative  procedure  was 
developed,  which  is  applied  to  the  nodes  in  the  fine  part  of  the  mesh: 
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First,  the  non-self  integrals  are  evaluated  at  all  nodes,  omitting  the  con¬ 
tributions  from  the  nearby  nodes.  Then  the  contributions  from  the  nearby 
nodes  are  added,  including  in  the  interpolation  procedure  the  newly  calcu¬ 
lated  field  values  at  the  current  time  points.  One  repetition  of  the  last 
step  gives  sufficiently  good  results.  The  use  of  a  local  fine  mesh  and 
iterative  solution  has  eliminated  this  instability  in  the  case  of  hard  and 
soft  targets.  At  time  of  writing,  this  procedure  has  not  been  included  in 
the  elastic  solution,  partly  in  order  to  minimize  complexity  at  this  ini¬ 
tial  stage  and  partly  because  other  instability  problems  were  felt  to  be 
more  severe. 

A  third  type  of  instability  can  be  identified,  named  the  resonance 
instability,  which  appears  to  be  caused  by  an  interaction  of  the  discrete 
formulation  of  space  and  time  points.  This  instability  is  not  in  phase 
on  the  surface  and  is  not  necessarily  oscillatory  (at  some  nodes  the  field 
amplitude  may  increase  exponentially) .  It  shows  up  in  the  far  field  cal¬ 
culation  as  an  exponentially  increasing  response.  This  problem  has  not 
been  solved,  but  it  appears  only  for  certain  unfortunate  choices  of  time 
step.  Slight  increases  or  decreases  in  the  time  step  can  eliminate  this 
instability. 

3.5.2  Surface  Derivative  Instabilities 

The  remaining  instability  problem  appears  only  in  the  general 
elastic  calculation.  It  has  the  nature  of  a  space-like  instability, 
indicating  that  it  is  due  to  the  inaccurate  determination  of  some  quantity, 
which  error  is  then  propagated  with  amplification  over  the  surface.  As 
can  be  seen  in  the  results  presented  in  the  next  section,  this  error 
appears  in  the  far  field  response  as  a  slow  oscillation  of  increasing  am¬ 
plitude.  This  problem  is  particularly  severe  because  its  onset  is  at  the 
beginning  of  time.  This  problem  has  not  been  solved.  The  likely  cause 

is  the  determination  of  the  surface  derivatives,  in  particular  of  the  term 
2 

V  u  .  It  should  be  noted  that  the  surface  derivates  are  taken,  not  of 


quantities  determined  directly  by  the  STIE,  but  of  u  which  is  obtained  by 

2-  2 

numerical  integration  of  9  u/9t  .  This  process  is  another  source  of  inac¬ 
curacy.  Proposals  for  solving  this  instability  problem  include  the  use  of 
the  iterative  solution  in  conjunction  with  a  finer  mesh,  and  the  formula¬ 
tion  of  the  problem  with  finite-element  techniques  which  would  avoid  the 
need  for  space  derivatives. 


3.6  THE  ANALYTIC  SOLUTION  FOR  A  SPHERE 


The  eigenvalue  solution  of  scattering  from  an  elastic  sphere  was 
implemented  in  order  to  test  the  STIE  solution  and  to  provide  another  means 
of  analyzing  the  measured  data.  This  solution  proceeds  by  writing  the  eigen¬ 
function  expansions  of  the  three  waves  and  solving  for  the  coefficients  by 
application  of  the  boundary  conditions.  The  result  is  a  very  complicated 
expression  for  the  coefficients,  which  will  not  be  duplicated  here,  but  is 
given  in  [1]  and  [2].  The  solution  (far  field  response)  is  written  in 
terms  of  frequency.  It  is  given  by  straightforward  evaluation  of  the  coef¬ 
ficients  and  of  the  Bessel  and  Hankel  functions  of  the  scattered  wave 
expansion. 

This  solution  was  also  used  to  examine  the  nature  of  the  response 
as  a  function  of  elastic  parameters. 

3.7  CONVOLUTION  PROCEDURE 

It  is  often  necessary  to  compare  responses  in  the  frequency 
domain,  or  to  compare  responses  with  different  inputs.  A  convolution  pro¬ 
cedure  accomplishes  this. 

Let  e(t)  be  the  input  signal  given  at  discrete  time  points,  and 
r(t)  the  discrete  response.  Then  the  complex  frequency  response  of  the 
target  is  given  by 

I  | 

H(oj)  =  R(co)/E(u))  ,  (53) 
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where  E(u))  =  r(t)j  ,  the  discrete  FFT,  and  where  the  division  represents 
the  discrete  (complex)  division  of  R  by  E  at  corresponding  values  of  to. 


The  time  domain  response  to  a  different  input  e^(t)  can  then  be 
calculated  by  forming 


r^(t)  =  F_1|h(oo)  •F(e'( t))j  , 

where  F  1  is  the  inverse  discrete  FFT.  Of  course,  H(to)  and  r^(t)  are  valid 
only  over  the  band  over  which  E(to)  is  defined.  The  magnitude  of  H  is 
plotted  in  the  results  to  follow. 

3.8  SCALING 


It  is  convenient  to  present  the  target  response  in  a  form  inde¬ 
pendent  of  scattering  geometry  and  medium  parameters.  Therefore  we  write: 

time  in  units  of  ct/a 

frequency  in  units  of  ka  =  2tt  fa/c 

scattered  pressure  in  units  r^p/a  per  unit  incident  pressure. 

Here 

c  is  the  medium  speed  of  sound, 

a  is  a  typical  target  dimension  (radius  for  the  sphere  and 
cylinder 

r^  is  the  distance  of  the  target  from  the  observation  point. 

The  frequency  response  H(ka),  as  calculated  in  Section  3.7,  is 
then  independent  of  incident  wave  form  and  amplitude,  observation  distance, 
target  size,  and  elastic  properties  of  the  medium.  The  only  quantities 
of  importance  then  are  target  shape  and  the  ratios  of  target  density  and 
sound  velocities  to  those  of  the  medium. 

[It  is  to  be  noted,  that  the  scaling  given  above  differs  from  the 
commonly  used  "reflectivity".  The  latter  is  defined  as  O/cross-sectional 
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area,  where  0  is  the  "cross-section"  or  cross-sectional  area  of  an  equiva¬ 
lent  sphere.  The  power  reflected  is  then  0/ 4irr2  or  Tra2/4irr2  =  (a/2r  )2 

O  o  o 

for  the  sphere.  That  is,  whereas  the  reflectivity  is  1  (at  high  frequency) 
for  a  sphere,  our  scaled  response  is  1/2.] 


•A 
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SECTION  4 

COMPUTATIONAL  RESULTS 


In  this  section  we  first  give  the  results  of  analysis  of  elastic 
sphere  responses  based  on  the  analytic  solution.  This  will  be  followed 
by  a  summary  of  results  of  STIE  calculations  for  the  elastic  and  fluid 
sphere  and  cylinder,  and  of  hard  and  soft  cylinders  near  or  on  a  bottom. 


4.1  ANALYTIC  RESULTS  FOR  A  SPHERE  -  EFFECTS  OF  ELASTIC  PARAMETERS 

In  Fig.  5  are  presented  the  responses  in  the  frequency  domain 
of  aluminum  and  brass  spheres.  The  elastic  constants  used  are 


(m/s) 

c2  (m/s) 

1  p  (g/cm3) 

Aluminum  6061-T6 

6180 

3090 

2.706 

Naval  Brass 

4400 

2100 

8.394 

Medium 

c3  =  1490 

- 

1.026 

Also  shown  in  Fig.  5  are  the  responses  of  hard  and  soft  spheres.  Figure  6 
illustrates  the  same  responses  in  the  time  domain,  for  a  smoothed  impulse 
with  a  =  1.  A  feature  shared  by  all  responses  are  the  resonances  at  in¬ 
tervals  of  6(ka)  =1.22,  which  are  characteristic  of  the  size  of  the  sphere. 
This  phenomenon  is  due  to  the  existence  of  creep  waves.  The  creep  wave 
path  length  is  longer  than  that  of  the  specular  return  by  (2+lT)a, 
hence  the  creep  wave  is  reinforcing  if  (2+7T)a  =  nX  =  ,  (for  n 

integer),  or  ka = 1.222  n.  Correspondingly,  in  the  time  domain  we  observe 
a  response  at  6(ct/a)  =  2  +  tt  after  the  specular  return. 

A  feature  distinctive  of  the  elastic  solid  response  is  the  exis¬ 
tence  of  deep  nulls  in  the  response.  We  observe  that  these  nulls  are 
spaced  at  about  6(ka)  =2.4  for  aluminum  and  at  6(ka)  =1.6  for  brass,  which 
is  approximately  in  the  ratio  of  the  shear  velocities  of  the  two  materials. 
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It  is  difficult  to  interpret  these  resonances.  They  correspond  approxi¬ 
mately  to  the  difference  in  time  traversal  from  the  point  at  0=45°  to 
-45°  around  the  back  of  the  sphere.  That  is,  a  creep  wave  at  the  shear 
velocity.  The  45°  point  is  where  from  physical  considerations  we  expect 
the  largest  coupling  into  the  shear  wave.  In  the  time  domain,  for  aluminum, 
the  above  6(ka)  =2.4  corresponds  to  S(ct/a)  =2.6,  and  indeed  a  response  is 
seen  at  that  interval  after  the  specular  return.  However,  the  correspond¬ 
ing  shear  creep  return  is  not  observed  for  brass  (6(ct/a)  =3.9).  For 
brass  in  the  time  domain  we  observe  a  wiggle  at  about  6(ct/a)  =2.9  after 
the  specular  return.  This  corresponds  approximately  to  the  delay  of  a 
shear  wave  starting  at  0=45°,  traversing  the  sphere  through  the  center, 
reflecting  from  -135°  and  coupling  out  at  -45°.  We  also  note  that  the 
initial,  specular,  return  is  not  much  lower  in  amplitude  for  the  hard  and 
these  elastic  targets.  That  is,  not  much  acoustic  energy  is  coupled  into 
the  interior  compressional  wave.  This  is  probably  the  reason  that  a  return 
due  the  interior  compressional  wave,  which  has  reflected  from  the  back 
side,  is  not  apparent  (it  would  show  up  at  6(ct/a)  =1  for  aluminum). 

To  illustrate  the  effect  of  the  elastic  parameters  on  the  re¬ 
sponse,  we  compare  the  aluminum  response  with  that  of  materials  with 
different  shear  velocities  c  but  the  same  p  and  c  (Figs.  7  and  8) .  We 

*  1  i- 

see  that  the  specular  part  of  the  return  is  changed  but  little,  as  might 
be  expected  since  there  is  little  coupling  into  a  shear  wave  near  the 
specular  point.  The  main  change  is  in  the  amplitude  of  the  wiggle  after 
t=0,  which  is  thus  clearly  associated  with  the  shear  wave.  In  Figs.  9 
and  10  is  shown  the  effect  of  a  change  in  P1  and  the  effect  of  a  change 
in  c  ,  keeping  the  other  two  parameters  the  same  as  for  aluminum. 

4.2  STIE  RESULTS  FOR  ELASTIC  AND  FLUID  TARGETS 

The  computational  results  for  a  sphere  using  the  STIE  formula¬ 
tion  is  given  in  Fig.  11  for  aluminum  and  in  Fig.  12  for  brass.  These 
are  the  smoothed  impulse  responses.  Also  shown  are  the  calculated  hard 
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(a)  p1  -  2.7 ,  c1  -  4430.  o  -  .0202 


ka 

(b)  p,  -  2.7.  c1  -  10330,  - 


k< 

(c)  p,  -  1.0,  c1  -  6200,  o  -  .3333 


FIG.  9  Effect  of  compressional  velocity  c^  and  density  p)  on  frequency 
response  of  sphere.  Medium:  c  =  1480,  p  =  1 ;  target:  C2  =  3100. 
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target  and  fluid  target  responses.  The  calculated  hard  target  response 
agrees  with  analytic  results  to  within  a  few  percent.  The  fluid  response 
is  for  a  target  with  the  same  density  and  compressional  wave  velocity, 
but  with  no  shear  wave.  No  theoretical  response  is  available  for  compari¬ 
son,  but  it  can  be  inferred  to  be  correct  on  the  basis  of  the  considera¬ 
tions  of  the  previous  section  and  comparison  with  the  results  in  ref.  [12] . 
The  results  for  the  elastic  target  are  compared  with  the  theoretical  result 
(dotted  lines)  and  illustrate  the  remaining  difficulty  with  the  elastic 
solution  as  discussed  in  Section  3.5.2. 

The  calculated  response  of  a  fluid  cylinder  at  axial  incidence 
is  given  in  Figs.  13  &  14.  This  is  for  a  cylinder  with  L/ a =2,  for  aluminum 
and  brass.  Superposed  on  the  plots  is  the  hard  target  response  (dotted 
line) .  As  is  the  case  for  the  sphere,  the  fluidity  only  changes  the 
amplitude  of  the  response. 


4.3  STIE  RESULTS  -  TARGETS  ON  BOTTOM 

We  present  here  the  smoothed  impulse  response  of  several  geome¬ 
tries:  the  sphere  and  right-circular  cylinder  at  several  depths  of  pene¬ 
tration  into  the  bottom  and  above  the  bottom.  These  calculations  were  made 
for  sound-hard  and  sound-soft  targets  on  sound-hard  or  -soft  bottoms.  In 
each  case  we  compare  the  response  with  that  of  a  target  in  free  water.  In 
all  cases  the  smoothed  impulse  width  constant  is  a^  =  1;  that  is,  the  pulse- 
width  (at  the  1.5%  points)  is  4/a^  =  4.  The  target  size  is  comparable  to 
this  pulsewidth  -  all  spheres  and  cylinders  considered  here  have  a  diameter 
of  2,  with  the  cylinder  length/diameter  =  3/1.  The  observation  point  in 
all  cases  is  the  same  as  the  source  point,  namely  r  =10  and  z  =  2  (height 
above  bottom) . 
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4.3.1  Spheres  Near  a  Flat  Bottom 

In  Fig.  15  we  plot  the  response  of  hard  spheres  at  several  ele¬ 
vations  with  respect  to  a  hard  bottom.  The  vertical  scale  differs  for 
each  plot  so  that  the  shapes  may  be  more  readily  compared.  Superposed  on 
each  plot  is  the  response  of  a  free  target  at  the  correct  relative  ampli¬ 
tude. 


We  identify  the  initial  positive  pulse  as  the  specular  return  and 
the  second  positive  pulse  as  the  creeping  wave.  Note  that  the  shapes  of 
all  the  responses  are  almost  identical  to  the  response  of  a  free  target  - 
only  the  response  from  a  sphere  3/4  embedded  in  the  bottom  (case  5)  has  a 
creeping-wave  to  specular-return  ratio  that  is  appreciably  different  from 
the  others.  The  amplitude  of  the  response  shows  a  simple  variation  with 
height  above  the  bottom:  increasing  as  the  target  approaches  the  bottom 
and  decreasing  again  as  the  target  penetrates  the  bottom. 

In  Fig.  16  the  response  of  soft  spheres  near  a  hard  bottom  is 

plotted.  Similar  to  the  hard  sphere,  the  effect  of  the  bottom  on  the  soft 

sphere  response  is  mainly  on  the  amplitude.  An  exception  to  this  observa¬ 
tion  is  the  positive  backswing  for  targets  which  are  entirely  above  the 
bottom  (case  1  of  Fig.  16) . 

For  a  soft  bottom,  the  source  and  its  image  have  opposite  signs. 
The  scattered  pressure  is  thus  a  difference  of  two  similar  signals,  the  net 
signal  being  due  to  the  differences  in  path  length  associated  with  each. 

We  thus  expect  the  return  from  a  target  near  a  soft  bottom  to  look  somewhat 
like  a  derivative  of  the  free  target  return.  This  is  borne  out  in  Figs. 

17  and  18  which  show  the  response  of  hard  and  soft  targets,  respectively, 

for  various  elevations  of  the  target  with  respect  to  the  soft  bottom.  In 
comparison  with  the  hard-bottom  case,  the  soft-bottom  returns  are  quite 
complicated.  The  signals  show  an  extra  peak  with  additional  structure  in 
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FIG.  1 5  Smoothed  impulse  response  at  various  elevations  -  hard  sphere,  hard  bottom. 
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FIG.  16  Smoothed  impulse  response  at  various  elevations  —  soft  sphere,  hard  bottom. 
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FIG.  17 
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the  tails.  The  ratio  of  the  peak  amplitudes  and  their  spacing  depend 
strongly  on  the  location  of  the  target.  Note  that  the  first  three  peaks 
(hard  target)  or  two  peaks  (soft  target)  remain  at  about  the  same  spacing 
for  all  elevations  of  the  target.  These  first  peaks  must  be  considered 
to  be  the  net  result  of  the  various  specular  paths  the  pressure  wave  can 
take  (via  target,  target-bottom,  bottom-target,  or  bottom-target-bottom) ; 
whereas  the  relative  timing  of  the  components  changes  but  little  in  these 
configurations,  that  change  has  a  large  effect  on  the  amplitude  since  the 
components  tend  to  cancel.  In  the  tail  of  the  response  we  see  a  signal 
that  increases  and  moves  rapidly  to  the  left  as  the  target  is  lowered  into 
the  bottom.  This  part  of  the  response  is  due  to  a  more  complicated  inter¬ 
action  of  bottom  and  target,  involving  creepage  along  the  target  or  bottom. 

4.3.2  Upright  Cylinder  Near  a  Flat  Bottom 

In  Fig.  19  is  shown  the  response  of  a  hard  right-circular  cylin¬ 
der  (RCC)  near  a  hard  bottom  at  elevations 

(1)  ^  of  its  length  above  the  bottom 

(2)  On  the  bottom 

(3)  ^  penetrating 

(4)  £  penetrating,  and 

(5)  3/4  penetrating. 

In  Fig.  20  is  shown  the  soft  RCC  with  the  same  configurations.  All  the 
cases  are  compared  with  the  response  from  a  free  target  at  broadside  inci¬ 
dence. 

Curve  1  (both  hard  and  soft  targets  -  target  above  the  bottom) 
shows  considerable  broadening  of  the  pulse  in  comparison  with  the  free  tar¬ 
get  response.  This  can  be  considered  due  to  the  large  separation  of  the 
target  and  its  image  (or  due  to  the  large  difference  in  the  various  paths 
of  specular  returns) .  Associated  with  the  broadening  is  a  reduction  in 
amplitude. 
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NOTE:  VERTICAL  SCALES  DIFFER-CHOSEN 
SUCH  THAT  TOTAL  SWINGS  ARE 
SAME  SIZE. 

DASHED  CURVES  ARE  FREE  TARGET  RESPONSE. 
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FIG.  19  Smoothed  impulse  response  at  various  elevations  -  hard  RCC  upright 
(a/L  =  1/6),  hard  bottom. 
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For  the  target  standing  on  the  bottom  (Curve  2),  the  pulse  nar¬ 
rows  and  increases  in  amplitude;  as  the  target  penetrates  the  bottom,  the 
response  decreases  in  amplitude  (due  to  the  decreased  effective  size  of 
the  target) .  It  should  be  noted  that  the  comparison  made  in  the  figures 
with  the  broadside  free  target  response  is  not  the  best  comparison,  since 
the  incidence  is  different  from  90°  in  all  the  targets  with  bottom  con¬ 
figurations.  Some  details  of  the  pulse  shapes  here  are  due  to  the  inci¬ 
dence  angle;  In  particular,  the  broadening  or  development  of  the  double 
bump  in  the  negative  swing  of  the  hard  target  response  (Fig.  19,  Curves 
1,2,3)  is  due  to  the  incidence  angle.  In  summary,  we  observe  that  as 
also  observed  for  the  sphere,  the  presence  of  a  hard  bottom  has  little 
effect  on  the  pulse  shape;  the  pulse  shape  being  determined  primarily  by 
the  effective  size  and  orientation  of  the  target  with  respect  to  the 
source.  The  small  doublet  wiggles  at  the  ends  of  the  responses  in  Fig.  20 
are  thought  to  be  caused  by  numerical  noise  in  the  solution.  They  become 
smaller  or  disappear  as  either  the  target  patches  are  made  smaller  or  the 
pulsewidth  is  increased. 

In  contrast  to  the  above  observation,  a  soft  bottom  causes  con¬ 
siderable  modification  of  the  return,  the  shape  of  the  pulse  being  strongly 
dependent  on  the  elevation  of  the  target  above  the  bottom.  See  Figs.  21 
and  22  for  the  hard  and  soft  right-circular  cylinders,  respectively. 

4.3.3  Cylinder  on  Its  Side  Near  a  Flat  Bottom 

The  next  configuration  evaluated  is  the  right-circular  cylinder 
lying  on  its  side  (i.e.,  with  its  axis  parallel  to  the  bottom).  In  Fig. 

23  are  shown  the  responses  for  a  "free"  cylinder  and  for  the  cylinder  above 
(center  at  1  diameter) ,  on  (tangent  to) ,  and  half-way  in  the  bottom.  At 
the  left  side  of  the  figure  the  target  is  viewed  in  line  with  the  axis, 
while  the  right  side  of  the  figure  shows  the  response  for  broadside  inci- 
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dence.  Again,  the  point  of  observation  is  the  same  as  the  source  point 
(r  = 10,  z  =  2)  .  The  free  target  response  is  shown  dashed,  superposed  on 
the  other  plots.  The  axially  viewed  response  consists  of  a  return  from 
one  end  followed  much  later  by  a  response  from  the  far  end,  while  the 
broadside  view  response  is  mainly  due  to  reflection  from  the  cylinder  side. 
We  see  that  the  amplitude  of  the  signal  increases  as  the  target  approaches 
the  bottom,  decreasing  again  as  it  penetrates  the  bottom.  An  interesting 
effect  is  the  near  disappearance  of  the  creep  return  for  the  target  tangent 
to  the  bottom,  viewed  broadside. 

Figure  24  illustrates  a  soft  target  with  a  hard  bottom  in  the 
same  configurations.  As  in  Fig.  20,  the  wiggles  at  the  ends  of  the  pulses 
are  due  to  numerical  noise  in  the  solution. 

Without  further  comment,  we  present  the  target  in  the  same  con¬ 
figurations  as  described  above  but  with  a  soft  bottom  in  Figs.  25  and  26. 
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FIG.  24  Smoothed  impulse  response  at  two  views  and  three  elevations  —  soft  RCC 
with  axis  horizontal  (a/L  =  1/6),  hard  bottom. 


FIG.  25  Smoothed  impulse  response  at  two  views  and  three  elevations  —  hard  RCC 
with  axis  horizontal  (a/L  =  1/6),  soft  bottom. 
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SECTION  5 

ANALYSIS  OF  MEASURED  DATA 


Measurements  were  performed  by  the  NCSL  on  a  variety  of  targets: 
sphere,  right  circular  cylinder,  aluminum,  brass,  several  angles  of  incidence. 

5.1  EXPERIMENTAL  PARAMETERS 

The  description  of  the  experimental  arrangement  as  reported  by 
NCSL  [13]  is  summarized  here. 

The  experimental  arrangement  is  shown  in  Fig.  27.  The  input 
signal  was  measured  at  the  target  location;  a  reverberation  signal  without 
target  was  measured  with  the  projector  and  receiver  at  their  normal  posi¬ 
tions  for  response  measurements.  The  projector  was  an  EDO  Western  6122, 
the  receiver  an  F-30.  The  alignment  of  the  target  axis  for  the  cylinder 
was  fixed  by  rotating  the  target  such  as  to  obtain  maximum  response  at 
axial  and  normal  incidence.  Each  measurement  was  repeated  100  times. 

The  speed  of  sound  in  water  was  determined  by  signal  delay 
measurement,  c  =  1491  m/s.  The  density  of  salt  water  was  assumed  to  be 

3  3 

1.026  g/cm  .  The  target  dimensions  were  a =  3  cm  for  the  spheres  and 
L/a  =  18  cm/3  cm  for  the  cylinders.  The  materials  were  Aluminum  6061-T6 
and  Naval  Brass  (Alloy  464) ,  for  which  the  elastic  parameters  are 


Aluminum 

Brass 

Y 

!  10 

15 

io6 

psi 

B 

!  10 

15.7 

io6 

psi 

° 

.33 

.341 

8.394  q/cm“ 


D 


2.706 


RESEARCH  CENTER 


(c) 

FIG.  27  Experimental  arrangement:  (a)  response  measurement,  (b)  input  measurement, 
and  (c)  reverberation  measurement. 
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The  desired  parameters  and  were  derived  from  the  relations  in  2.1.3. 
For  aluminum  this  gave  reasonable  results.  For  brass,  the  constants  on  the 
preceding  page  were  made  self-consistent  by  using  0  =  0.341  rather  than  the 
value  of  .331  quoted  in  materials  handbooks.  Anticipating  the  results  of 
the  analysis,  it  was  found  that  good  agreement  with  measurement  was  obtained 
if  c2  is  taken  as  2%  smaller  than  from  direct  use  of  Y  and  B  in  the  equa¬ 
tions  of  Section  2.1.3.  (It  was  assumed  that  this  adjustment  was  justified 
in  the  light  of  the  inaccuracy  in  the  preceding  parameters.)  Accordingly, 
the  parameters  used  for  analysis  were 
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A  sample  of  raw  data  (aluminum)  sphere  is  given  in  Fig.  28. 

5 . 2  DATA  REDUCTION 

The  data  was  put  into  the  standard  form  discussed  in  Section 
3.8,  by  scaling  time  to  ct/a,  frequency  to  ka,  and  the  response  amplitude 
by  r^/a.  The  amplitudes  of  the  response  were  adjusted  by  the  quoted  re¬ 
ceiver  gain.  (Also,  for  convenience  in  plotting,  the  amplitudes  of  all 
signals  were  scaled  by  an  equal  factor,  F,  which  disappears  when  the  fre¬ 
quency  response  is  computed.) 

The  100  individual  pings  of  the  measurements  were  found  to  be 
extremely  similar  in  phase,  and  in  amplitude  differed  only  with  a  standard 
deviation  of  2%  for  the  responses  and  10%  for  the  reverberation.  The  100 
pings  of  each  measurement  were  simply  averaged,  discarding  only  those  (few) 
pings  which  showed  spurious  behavior. 
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Figure  28  illustrates  an  example  of  scaled  but  unprocessed  meas¬ 
urements  (for  the  brass  sphere  at  0°).  The  incident  pressure  is  taken  as 
the  voltage  at  the  receiver,  vgl  (since  the  receiver  was  placed  at  the 

target  position  for  this  measurement) ,  scaled  by  the  receiver  gain,  G  , 

e 

and  an  arbitrary  factor,  F.  That  is: 

-(G  /20) 

inn  e 

e (ct/a)  =  p  (ct/a)  =  F  •  10  v  (ct/a) . 

e 

The  measured  response,  which  consists  of  the  actual  response,  r,  plus 
noise,  n,  is  scaled  similarly  and  also  by  the  factor  rQ/a.  That  is: 


r  -(G  /20) 

r+n(ct/a)  =  —  •  F  •  10  m  v  (ct/a). 

6l  in 


The  scaled  reverberation  noise  signal  is  represented  as 


r  -(G  /20) 

n(ct/a)  =  —  *  F  •  10  n  v  (ct/a)  . 
a  n 


These  same  signals,  transformed  to  the  frequency  domain,  are  illustrated 
in  Fig.  29. 


The  scaled  response  with  noise  subtracted  is  then 

r  /  -(G /20)  -(G  /20) 

r(ct/a)  =  — —  *  F  •  I  10  v  (ct/a)  -  10  n  v  (ct/a) 

sl  \  m  n 

It  was  expected  that  direct  subtraction  of  the  reverberation  signal  would 
yield  a  clean  response  signal.  This  was  not  the  case.  Some  reverberation 
noise  remained  in  the  precurser  and  after  signal  end.  The  result  of  noise 
subtraction  is  shown  in  the  first  curve  of  Fig.  30.  (Arbitrary  adjustment 
of  reverberation  amplitude  did  not  affect  the  overall  result.)  A  possible 
cause  of  this  problem  is  environmental  changes. 
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FIG.  29 
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Typical  response,  input,  reverberation  in  frequency  domain  (brass  sphere  0°) 
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The  effect  of  the  remaining  noise  is  to  introduce  spikyness  in 
the  frequency  response.  To  reduce  this  effect,  smoothing  was  introduced 
by  the  application  of  a  cosine  bell  window,  w,  to  the  signal,  which  elimi¬ 
nated  the  initial  and  final  reverberation.  Fig.  30(b)  shows  the  signal 
after  reverberation  subtraction  and  application  of  the  window,  and  also 
shows  the  window  used.  Figure  31  illustrates  these  results  of  signal 
processing  in  the  frequency  domain.  The  (complex)  frequency  response  of 
the  target  was  then  found:  H(ka)  =[k(ka)  *  W  (ka)]/E  (ka)  .  The  result  of 
this  process  is  shown  in  Fig.  32  for  the  aluminum  sphere  and  in  Fig.  33 
for  the  naval  brass  sphere.  The  results  are  shown  for  the  two  aspects 
(0°  and  45°)  at  which  the  spheres  were  measured.  The  plotted  quantity  is 
the  magnitude  of  H(ka).  Note  that  the  input  signal  occupied  a  band  of 
ka=1.5  to  12,  so  that  H(ka)  outside  these  limits  has  no  meaning. 

5.3  COMPARISON  WITH  CALCULATIONS 

5.3.1  Spheres 

A  slight  directionality  can  be  observed  in  comparing  the  0°  and 
45°  aspects,  indicating  either  a  slight  anisotropy  of  the  material,  or  an 
error  due  to  changes  in  the  reverberation.  The  latter  is  likely,  since 
the  reverberation  spectrum  was  greatest  at  those  frequencies  where  the 
differences  are  observed. 

Figures  32  and  33  also  show  the  theoretical  responses  superposed 
as  dotted  curves.  The  results  are  in  close  agreement  with  respect  to  both 
amplitude  and  alignment  of  the  major  nulls.  (The  value  of  c for  brass 
was  adjusted  by  2%  as  noted  in  Section  5.1.)  The  faster  small  wiggles 
with  null  spacing  at  6 (ka)  =1.22  does  not  agree  as  well,  but  is  clearly 
present  in  the  response.  The  small  and  deep  wiggles  can  thus  be  used  to 
estimate  the  size  and  composition  of  the  target,  independently  of  size 
determination  based  on  response  amplitude.  This  is  perhaps  even  more  appa¬ 
rent  if  the  signals  are  observed  in  the  time  domain  (see  Section  4.1) 
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(a) 

45°  ASPECT 

/ 

A 


THEORETICAL 


which  yields  directly  circumferential  size  from  the  creep  recurn  and  target 
material  from  the  wiggle  in  the  response  near  t=0.  Figure  34  illustrates 
the  theoretical  responses  in  the  time  domain  for  the  two  spheres  calculated 
by  convolving  the  input,  e,  with  the  theoretical  H(ka): 


r'(ct/a)  =  F  1|F^e(ct/a))  •  H(ka) 

The  agreement  with  measurements  is  seen  to  be  close,  differing  somewhat  in 
amplitude,  but  not  in  zero  crossings. 


5.3.2  Cylinders 

The  frequency  responses  of  the  cylinders  at  0°  and  90°  are  shown 
in  Figures  35  and  36  for  aluminum  and  brass  respectively.  Superposed  on 
these  signals  is  the  high  frequency  asymptotic  response.  For  09  (end-on) 
aspect  the  latter  is  given  by 

H (ka)  =  0.5  (ka) 

for  90°  (broadside)  aspect  by 


H  (ka) 


^  k  (ha)1/2  . 

y?  a 


It  can  be  seen  that  the  high  frequency  approximation  is  not  a 
good  model  for  the  elastic  cylinders  in  this  frequency  range. 
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(b)  BROADSIDE  ASPECT 


FIG.  35  Measured  response  of  cylinder  (aluminum). 
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(a)  END-ON  ASPECT 


(b)  BROADSIDE  ASPECT  7s-tas 

FIG.  36  Measured  response  of  cylinder  (brass). 
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SECTION  6 

FORMULATION  FOR  HULLED  BODIES 


We  consider  now  the  problem  of  a  fluid-filled  elastic  shell,  sub¬ 
merged  in  a  fluid.  Following  the  formulation  of  Section  2  and  writing  u^ 
for  the  interior  fluid,  u2  for  the  shell,  and  u3  for  the  exterior  fluid, 
we  thus  have  the  wave  equations: 


in  the  fluids,  and 


(V  •  u.)  =  c2  V2(V  •  un) 
3t  1  1  1 


(54) 


3t 


\  (V  •  u3)  =  c2  V2(V  •  u3) 


(55) 


^2  (V  •  u2)  =  C21  V2(V  *  G2} 


3t 


(56) 


~~2  (V*V  =  C22  V2(VX  V 

dt 


(57) 


in  the  solid  shell,  where  c^  and  c ^  are  the  shell  compressional  and 
shear  velocities  respectively. 


These  equations  can  be  solved  in  the  fashion  presented  in  Sec¬ 
tion  2.2.2  by  writing  the  space-time  integral  equations  for  the  boundary 
points.  The  result  is  an  integral  equation  for  (V  •  u^)  integrated  over 
the  interior  boundary,  an  integral  equation  for  (V  •  u^)  integrated  over 
the  exterior  boundary,  and  integral  equations  for  (V  •  u2)  and  (Vx^) 
integrated  over  the  sum  of  the  interior  and  exterior  boundaries.  These 
STIE's  are  then  solved  simultaneously  by  using  the  boundary  conditions  of 
Section  2.1.4*.  one  set  of  boundary  conditions  for  each  surface. 
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Such  a  "full-scale"  solution  thus  requires  the  evaluation  of  4 
integral  equations,  the  use  of  boundary  relations  on  both  interior  and 
exterior  shell  surfaces,  and  the  storage  of  11  field  variables  (2  more 
than  for  the  case  of  the  elastic  solid  in  a  fluid) .  The  solution  process 
is  thus  a  straightforward  extension  of  the  solution  for  the  elastic  solid 
problem. 


Alternatively,  the  problem  can  be  considerably  simplified  if 
the  shell  can  be  treated  as  a  thin  membrane.  By  "thin"  is  meant  that  there 
is  no  variation  in  u  in  the  normal  direction  (i.e. ,  across  the  membrane) 
and  that  the  only  stresses  in  the  shell  are  tangential  or  surface  tensions. 
The  first  boundary  condition  (see  Section  2.1.4)  then  becomes 


u  =  u  =  u  (continuity  of  normal  displacement) .  (58) 

nl  n2  n3 

The  second  boundary  condition  (continuity  of  normal  stress)  holds  indepen¬ 
dently  on  the  interior  and  exterior  surface  of  the  membrane,  but  not  neces¬ 
sarily  from  interior  fluid  to  exterior  fluid,  as  the  thickness  of  the  mem¬ 
brane  is  reduced  to  zero.  The  membrane  permits  a  discontinuity  of  normal 
stress  (i.e.,  pressure)  to  exist  between  the  fluids. 


The  forced  wave  equation  for  the  normal  displacement  of  a  thin 
membrane  is  [8] 


Here , 


a2 
d  u 


(P2d) 


n  „  n2 

r T  V  u 


3tz 


=  P1  '  P3 


(59) 


d  is  the  thickness  of  the  membrane, 

(p^d)  is  the  density/unit  area, 

T  is  the  surface  tension  (force/unit  length) , 
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Vg  is  the  surface  Laplacian, 

P-^  “  P3  is  the  discontinuity  in  excess  pressure, 
n  is  directed  into  medium  3  (or  outward) . 

The  surface  tension,  T,  is  in  general  a  function  of  the  surface 
coordinates.  For  a  membrane  that  is  free  to  move  in  the  fluid,  T  is  pro¬ 
portional  to  the  difference  in  ambient  pressures  between  media  1  and  3  and 
to  the  equilibrium  surface  curvature. 

Writing  our  equations  in  terms  of  excess  pressure  p,  rather  than 
(V  •  u)  (see  Section  2.1.5),  we  can  thus  summarize: 


Wave  equations: 


2  „2 

P1  =  C1  V  P1  ' 


••  2  _2 
p  =  c  V  p  , 
*3  3  p3  ' 


which  have  the  STIE  solutions  on  the  boundary  (see  Section  2.3.1): 


P1(r,t) 


_  _1_  f  (  1  _^1 
2tt  J  I  R  dn' 


,  2,  r  >  I  ,  _  > 

(n  •  R)  L1  px  j  dS  , 


P3 (r ,t)  =  2p 


-i  f  |i 

2tt  J  Jr 


3--  (r-R)  L^jdS'  . 


These  are  solved  simultaneously,  using  the  condition  (59)  and  the  relations 
(see  Section  2.1.5) 
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We  see  that  the  solution 
and  the  integration  of  ii 


PlUn 


P3Un 


(64) 


also  requires  the  formation  of  V2u  on  the  surface 

s  n 

in  time  to  give  u  . 

n 
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SECTION  7 
SUMMARY 

In  this  report  we  have  formulated  in  detail  the  general  solution 
to  time  domain  acoustic  scattering  from  elastic  targets,  both  solid  and 
thin-shelled.  The  formulation  was  implemented  for  solids  with  rotational 
symmetry  at  arbitrary  angle  of  incidence.  The  results  for  the  hard,  soft, 
and  fluid  approximations  were  very  good  for  both  smooth  (spherical)  bodies 
and  bodies  with  edges  (right  circular  cylinder) .  The  full  elastic  solution 
still  suffers  from  numerical  instabilities  believed  to  be  caused  by  inac¬ 
curacies  in  the  formation  of  surface  derivatives.  Also,  it  is  believed 
that  special  treatment  is  required  at  sharp  edges  for  the  full  elastic 


The  responses  of  spheres  and  right  circular  cylinders  on  or  near 
a  bottom  were  calculated  for  the  cases  of  hard  and  soft  targets  and  bottoms 
A  method  of  images  was  employed. 

The  responses  of  elastic  spheres  were  calculated  according  to  the 
classical  solution.  Using  this  solution,  the  effect  of  elastic  parameters 
on  the  response  was  analyzed. 


The  results  of  measurements  taken  at  NCSL  were  compared  with 
calculations.  It  was  found  that  the  hard,  fluid,  elastic  target  models, 
in  that  order,  result  in  increasingly  good  approximations  to  the  measured 
responses  for  incident  waves  in  the  resonance  region. 
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APPENDIX  A 

REFORMULATION  OF  THE  SPACE-TIME  INTEGRAL  EQUATIONS, 
USING  INTEGRATION  BY  PARTS 


An  alternate  approach  to  solving  the  STIE's  using  integration  by 
parts  is  described.  This  method  requires  the  storage  of  only  4  or  5  sur¬ 
face  variables  (rather  than  9) ,  at  the  expense  of  greater  complexity.  In¬ 
stability  difficulties  led  to  abandonment  of  this  method  in  favor  of  the 
less  complex  approach  in  the  body  of  the  report. 

We  start  here  with  equations  (24) -(26)  and  divide  the  integral 
into  a  self  and  non-self  part.  For  equation  (24)  we  had  obtained 


~2 
9  u 


(1  -  e) (V  •  u3)  +  y 


n 


c39t 


2  **  2(V-u3) 


-  .  inc  1  r  (  1 

~  J 


r\  2  > 

d  u 


2TT 


NS 


RcT  9t 


(R  •  n")  L  (V  •  u')[  dS' 

j  3  ) 


(A-l) 


We  plan  to  store  the  quantities  appearing  under  the  integral  u  ,  V  •  u 

n  3 

For  equation  (25) ,  we  proceed  similarly,  using 


d_ 

9n 


(V 


u) 


(A-2) 


(derived  from  B.C.  4  of  Section  2) ,  and  using  B.C.  2  to  eliminate 
V  •  u  from  the  equation.  We  then  obtain  for  the  self-integral: 


(A- 3) 
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(V  •  u) 


(A-4) 


The  non-self  integral  (which  has  the  same  argument  as  in  equation  (A-3)), 

contains  in  addition  to  u  and  V  •  u  ,  the  surface  derivatives  of  u.  We 

n  3 

wish  to  eliminate  these,  using  integration  by  parts,  which  will  transform 
the  surface  derivatives  into  time  derivatives.  We  use  the  integration  by 
parts  formula  for  an  open  surface: 


J "  <P  (Vg  *  A)  dS  =  -  J  nc  •  <pA  dJL  -  J *  A  •  V<J>  dS 


(A— 5) 


Here,  C  is  the  contour  bounding  the  surface,  S,  of  integration  and  n  lies 

c 

in  S,  normal  to  C  pointing  into  S  (i.e.,  pointing  out  of  self-patch).  For 
example,  identifying  A  as  u  ,  <|>  as  (R  •  n  )L  (where  the  subscript  S  denotes 
surface  components) ,  the  term  in  the  non-self  integral: 

(4-)  W='lds' 

NS  C1  ' 


becomes 
2 


2  q  2  * 

T  *  H  K2)  /  ’  S'>  VSc  •  \’d*'  ♦  Jn  (-2-)  /  «  •  «'>  VS  • 

1  c 


dS' 


NS 


(A-6) 


where 


0 i  is  the 


operator  (  +  — - - —  \ 

\  R3  Re?  9t  rc?  9t2  ) 


(A-7) 
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This  is  obtained  from  the  identity 


V>  =  =  (R.n 'IRO^^  .  (A-f 

The  contour  integral  about  the  self-patch  can  be  evaluated,  using  a  Taylor 
series  expansion,  and  results  in  terms  similar  to  those  in  equation  (A-4) , 
and  additional  terms  involving  the  difference  in  principal  curvatures. 

Similarly,  the  non-self  integral  of  f  V2  u'  ds'  can  be  evalua 

J  s  n 
NS 

ted  by  applying  the  integration  by  parts  formula  twice.  The  final  result 
for  the  transformed  V  •  u  equation  is 

o  c2  2  a2 


c~'J  c-at' 


Isl  + 


2  '  2^2 
C1  Cl9t 


+  (  — rr~)  (R  •  n ')  0  (R  •  u")j 

'  C;  ’  1  8  I 


where 
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Only  the  quantities  u  and  (V  •  u^)  appear  under  the  integral  sign,  which 

(four)  quantities  need  to  be  stored.  The  terms  in  I  ,  need  to  be  evalua- 

sl 

ted  by  numerical  surface  differentiation  of  u.  Equations  (A-9)  and  (A-l) 

—  2  2 

can  then  be  solved  simultaneously  for  (V  •  u,)  and  8  u  /3t  .  Furthermore, 

3  n 

using  the  representation 


3  u  u  (t  +  At)  -  2  u  (t)  +  u  (t  -  At) 

_ n  ...  _  _n _ n _ n _ 

2  2 
3t  (At) 


(A- 10) 


u^ft)  can  be  evaluated  for  the  next  time  point. 


In  a  similar  manner,  the  (V  x u)  equation  can  be  transformed  to 
2  _  2 

yield  3  (nxu)/3t  and,  by  use  of  (44),  u  and  u  at  the  next  time  point. 

x  y 

The  normal  component  of  the  (V  x  u)  equation  yields  the  normal  component  of 

(V  x u) .  Five  quantities  appear  under  the  integral  sign:  (V  •  u^) ,  u,  and 

(V  x  u)  . 
n 

The  equation,  which  has  a  lengthy  derivation,  is  written  here 
without  proof: 

<l*e><v,S>  S-X-A^a,;)  =Is;,*l8s  (a-11) 

c29t 


where 


?s2  =  -(1+f)  (VxG)s  -  fix  V(V  *  S’  +  ^  ‘  f  Y 


K„  -  K„v  ,  ,3u  3u 


-iYfix7«7s.5,  -  y( 


3u  \  3u  3u 

_ y  \  ,  n  ,  ~  i 

+  ^r)+  x  i7  +  y  ii 


3u  v 

IT  • 
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The  solution  procedure  can  thus  be  summarized: 

at  time  t, 

(1)  Evaluate  the  non-self  integrals  for  the  three  equations 

(2  scalar  and  one  vector  equation) ,  using  the  retarded 

values  of  the  five  stored  quantities  (V  •  u„) ,  u,  (V  x u)  . 

3  n 

(2)  Evaluate  the  self-terms  involving  surface  derivatives  of 
u(t) ,  where  u(t)  had  been  established  at  the  previous  time 
iteration. 

(3)  Solve  (A-9)  and  (A-l)  for  V  *u,(t)  and  u  (t  +  At) 

J  n 

(4)  Use  7  •  u,(t)  to  evaluate  the  remaining  term  in  I  .  of  (A-ll) 

3  s2 

(5)  Solve  (A-ll)  for  (Vxu)  (t)  and  u  (t+At),  u  (t  +  At). 

n  x  y 

The  process  is  now  repeated  for  the  next  time  point. 


The  quantity  (7  x  u)  can  also  be  found  by  surface  differentiation 
n 

of  u^  and  u^,  and  hence  can  be  eliminated  entirely  from  all  the  equations. 

It  is  thus  possible  to  solve  the  problem  using  only  four  stored  quantities. 

(The  only  reason  (7  x  u)  is  calculated  and  stored  in  this  formulation  is 

n 

that  it  is  anticipated  that  this  quantity  may  give  rise  to  difficulties 
in  the  solution,  since  it  is  overdetermined.)  We  might  expect  from  physi¬ 
cal  considerations  that  four  quantities  is  all  that  is  required:  Inside 
the  solid  the  vector  field  u  (which  satisfies  a  wave  equation)  is  specified 
by  6  quantities  on  the  surface  (for  instance,  u  and  its  normal  derivative); 
outside,  in  the  fluid  medium,  the  scalar  field  (7  •  u^)  or  Pe  is  specified 
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by  2  quantities  on  the  surface.  Now,  the  boundary  conditions  impose  4 
relations  among  these  quantities,  leaving  a  net  of  4  independent  quantities 
It  is  seen  above  that  the  equations  involving  only  4  quantities  are  con¬ 
siderably  more  complicated  than  the  more  natural  formulation  involving  10 
quantities. 


